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^ ■ Abstract 

We point out several subtleties arising in brane-world scenarios of cosmological constant 
cancellation. We show that solutions with curvature singularities are inconsistent, unless the 
contribution to the effective four-dimentional cosmological constant of the physics that resolves 
the singularities is fine-tuned. This holds for both flat and curved branes. Irrespective of this 
problem, we then study an isolated class of flat solutions in models where a bulk scalar field 
with a vanishing potential couples to a 3-brane. We give an example where the introduction of 
a bulk scalar potential results in a nonzero cosmological constant. Finally we comment on the 
stability of classical solutions of the brane system with respect to quantum corrections. 



1 Introduction 



With the discovery of string dualities it became clear that extended objects (branes) are built 
into any string theory. In the form of D-branes [0 they first served as tools in addressing concep- 
tual questions of string theory (e.g. dualities, black hole entropy, holography). An interesting 
property of D-branes is that they provide a natural picture for having some fields (typically 
gauge fields) being confined to a hypersurface in space whereas others (typically gravitational 
fields) can propagate in all directions. This property makes contact to an early suggestion 
|| [| that the matter in our universe is confined to live on a 3+1 dimensional hypersurface 
of a higher dimensional universe. In the simplest case the geometry is such that the higher 
dimensional space is a direct product of our 3+1 dimensional universe with some additional 
internal space. This picture was generalized to warped compactifications [|J, where the 3+1 
dimensional metric depends on the position of the brane in the higher dimensional space. The 
authors of M constructed a model where this dependence is exponential. They were led to the 
conjecture that a natural explanation for the large hierarchy between the electroweak scale and 
the Planck scale was found. But even without solving the hierarchy problem warped compacti- 
fications open up interesting possibilities, like localizing gravity on the brane with an infinitely 
large extra dimension ||. 

Apart from the hierarchy problem it has been also tried to solve the cosmological constant 
problem^] within the brane world scenario ||, pLOf] - The basic new ingredient is that there is a 
scalar living in the bulk with coupling to the brane. Naively one then hopes that the well-known 
fine-tuning of the parameters of the theory required in order to ensure a vanishing cosmological 
constant is replaced by the adjustement of the zero mode of this additional scalar. However, 
these models typically have singularities being located within a finite distance from the brane. 



Problems arising due to the singularities have for example been pointed out in Jy], (this and 



other aspects of these models have been also the subject of Refs. [|T3|, [14], [15], [16], [T7|, [L8|, [LJ|). 



In our note |12[] we observed that the singularities can be given a physical interpretation in the 
presence of additional source terms. Then the curvature singularity just reflects the fact that 
we are considering the limit in which the sources (or branes) at the singularities have no finite 
extension into the transverse direction (this is similar to black hole singularities appearing due 
to the point like nature of the source). With these additional sources, however, a fine-tuning 
is needed to obtain a vanishing cosmological constant for the effective four dimensional theory. 
Irrespective of this problem, it has been stressed in p, |2(J that for a very particular choice of 
parameters (namely vanishing bulk scalar potential and particular form and value of the scalar- 
brane coupling) there exist isolated flat solutions, i.e. there seems to be no smooth deformation 
connecting zero cosmological constant to a non-vanishing value. 

In the present paper we are going to elaborate on brane world models as a framework to find 
solutions of the cosmological constant problem. In particular we want to clarify the following 
points: 

• the vanishing of the cosmological constant requires a consistency condition for the brane- 
world setup which is especially crucial in the presence of singularities (this has been 
sometimes overlooked in the literature); 



1 For reviews of the cosmological constant problem see @, 0, §• 
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• known mechanisms to fullfill this consitency condition require a fine tuning of parameters 
of the model, comparable to the usual fine tuning of the cosmo logical constant. 

The so-called self tuning solution of the cosmological constant problem is therefore at best a 
scenario where this problem has been rephrased. The real solution of the problem would have 
to explain the fine tuning that is necessary for the consistency condition mention above. 

The paper is organized as follows. In the next section we specify the general setup and 
recall consistency conditions on warped compactifications. In section three we generalize the 
points made in We show in various examples (flat and curved branes) that one can make 
solutions with singularities consistent by adding additional branes to the model. To achieve a 
definite value for the cosmological constant (e.g. zero) a fine-tuning is needed. Section four is 
devoted to a careful investigation of the isolated flat solution mentioned above, for which we 
find that introducing a non trivial bulk potential for the scalar field necessarily results in a non 
zero cosmological constant. In the fifth section we discuss the relation of the classical solutions 
of the brane-bulk system to the physical quantum world. Finally the last section is devoted to 
a summary and outlook. 

2 General Setup 

To study thoroughly the structure of vacuum solutions in the brane worlds it is useful to 
consider an ansatz corresponding to maximally symmetric 4d foliations, which beyond the 
Minkowski space include de Sitter and anti-de Sitter 4d space-times. The simplest nontrivial 
setup corresponds to just a single transverse coordinate. The line element is 



with g^ v = g^ v + where the background metric corresponds to one of the 4d maximally sym- 
metric spaces: g^ = diag (— 1, +1, +1, +1) (Minkowski metric); diag (— 1, e 2v/ ^*, e 2v ^', e 2v ^-^ 

((IS4 background with A > 0); diag ( - e 2 V^ x \ e 2yf-Kx 3 ^ e 2y/^Kx 3 ^ ^ \ (^yg 4 background, A < 



0). With we denote a small fluctuation around the background metric. The sources for 
such configurations are assumed to be located on a number of four dimensional branes, one 
of which represents the observable gauge sector. It is important to note that although the 
observable gauge interactions are strictly confined to the 3-brane, the gravity and moduli fields 
permeate the whole space, effectively connecting the walls in a nontrivial way. This implies 
that every particle localized on the wall feels sources of gravitational forces which are located 
all over the bulk. In some cases the influence of the remote sources will be suppressed, like 
in the case of the exponentially falling off graviton wave function in ref. 0, which would 
effectively restrict the relevant gravitational sources to the thin layer around the brane, but 
sometimes the suppression would be so mild, that the influence of the whole bulk contribution 
will be highly relevant. Thus in any case, even though the gauge forces are restricted to the 
branes, the gravitational sector has to be completely integrated out when going to the effective 
four dimensional theory. In particular, this implies that when one computes the effective four 
dimensional energy density, or four dimensional vacuum pressure, one has to integrate over the 
whole causally accessible portion of the transverse space. 




(1) 
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Let us analyze the relation between 4d and 5d physics in some detail. The 5d action we 
consider is (we follow the conventions of Ref. [Kj for the normalization of the Einstein term 
and of the scalar kinetic term) 

S 5 = J d 5 xV^G (r - ~ (d<j)) 2 - V (0)) + J m d 5 xV=9 (-fi (0) 5 (x 5 - xf)) (2) 

where V(4>) is the bulk scalar potential, fi(<fi) are brane contributions to the Lagrangian, the 
index i counts the branes, and g^ u is the induced metric on the branes. The corresponding 
Einstein equations are 

1 



V^G[E MN --T MN J = , (3) 
where E MN = Rmn — \ GmnR is the Einstein tensor and 

T MN = I d M <t>d N <f) - ^ (d<p) 2 G MN - V(<p)G MN - f t 5 (x b -xf) g^ M 5 v N (4) 
is the energy-momentum tensor. The dilaton equation of motion is 

_ ^-^G - (x 5 - x 5 A + ~d M (V^GG MN d N <f>) = 0, (5) 



where M, N = 1, . . . , 5 and /i, v = 1, . . . , 4. Although we do not consider fluctuations of G 55 
that would correspond in four dimensions to the modulus associated with the length of the 
fifth dimension and put here G55 = 1, we shall comment on the stability of the fifth dimension 
in the case of static solutions at the end of the paper. With the choice g^ u = T]^, the ansatz 
([!]) naively implies that the 4d cosmological constant should vanish. As explained in |T2|] this 



is the case if the ansatz solves Einstein equations globally, i.e. everywhere in 5d space-time. 
If however there happens to be singularities, for instance at points where = 0, then 

one needs to 'repair' the system at such points in order to achieve truly flat 4d background. 
This issue shall be discussed later on. In the presence of the general sources the warp factor, 
hence also the curvature scalar, as well as the bulk scalar do assume some nontrivial and 
^/-dependent vacuum configuration. 

To identify 4d gravity and the strenght of its coupling to matter one needs to split the 
action (||]) into vacuum part, and fluctuations around it. The vacuum part is read off from the 
5d action upon substituting solutions of the equations of motion for the warp factor and for 
the scalar. This vacuum Lagrangian is 

(Lb) = Jdye< A ^)^-Vm) + \{h)5{x"-xf ) ) ) (6) 

where (...) indicates that the classical background value of the corresponding field is taken. It 
is immediate to see that the lowest terms in of the expansion of the action (pD around the 
vacuum are 

S 5 = d 5 xJ^M 3 e 2A ^R+... (7) 



where R denotes the curvature built out of g^ u . From (0) one can identify the 4d Planck scale 
Mp t = M 3 J dye 2<A<yX - )> . To obtain the effective four dimensional action S4 for gravity one 



3 



needs to integrate over the fifth coordinate. The model is consistent if {g^J) from our ansatz 
([I]) minimizes 5*4. Thus the effective four dimensional action reads 



S 4 = M 2 Pl J d 4 x^g (R-\ 



where the cosmological constant A = 6A. Finally we arrive at a consistency condition by the 
requirement that the on-shell values of S 5 and S4 should be equal[| 

(L 5 > = M Pl ((R) - A) = 6AM Pl . (9) 

A breakdown of this condition would signal that Einstein equations are not satisfied everywhere 
in space-time. It is interesting to note that Eq. (§) can be rewritten 

dx 5 e 4<A> (T ° + T 5 5 ) = 6AM Pl . (10) 

In the case of a Poincare-invariant background, the vanishing of the 4d cosmological constant 
is therefore equivalent to a constraint on the 5d energy-momentum tensor, namely the integral 
on the left-hand side of ([H]) should vanish. This constraint is a variant of the condition 
Jdx 5 e <A> (T$ - ±Tf) = derived in Ref. |2T|, and it has the same origin: the combination 
e 4<A> (To + T|) is constrained by Einstein equations to be a total derivative, therefore its 
integral over a compact interval should vanish. However, while / dx 5 e <A> {TQ — |T|) = also 
holds for a dS± or AdS± background, the condition / dx 5 e 4<A> (T® + T|) = is more specific 
for Poincare invariance. To see how Eq. ([IT]) (with A = 0) relates the vanishing of the 4d 
cosmological constant to the consistency of the 5d Einstein equations, let us note that in the 
Randall-Sundrum model g r c J dx 5 e 4<A> (T ° + T|) = (e" 4 ^ rc - 1) A/A; - e- 47Tkr ^V obs - V hU ; 
therefore, the above constraint is automatically satisfied once the fine-tuning needed to achieve 
a Poincare-invariant solution in four dimensions (namely Vhid — —V bs = A/k) is imposed. A 
similar cancellation is expected to appear in the approach put forward in Ref. 0, |22], p3| . 



3 Vacuum configurations with singularities and fine-tuning 

In this section, we consider various vacuum configurations leading to singular warped compact- 
ifications, and show that in this case Einstein equations are not globally satisfied (condition 
([Top is not fulfilled), leading to an inconsistency of the solution. This implies in particular that 
the recently proposed "self-tuning" mechanism of the cosmological constant || [lCj does not 
work as it stands, and hides a fine-tuning at the singularity. 

Let us first recall what is generally understood under "self-tuning of the cosmological con- 
stant" in brane-world models. It is the property that solutions of the five-dimensional Einstein 
equations preserving Poincare invariance on the brane can be found for a wide range of values 
of the brane tension, i.e. independently of the brane contribution to the effective cosmological 



constant. This behaviour, first identified in Ref. M and [10], has been shown to be generic in 



2 Note that for A = this is just the condition of vanishing vacuum energy imposed in (lj 
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models where a single brane is coupled to gravity and a scalar field flt^| . Indeed, plugging the 
ansatz 



ds 2 = e 2A{x5) r ]flu dx^dx v + (dx 5 f (11) 

into the equations of motion (|3]) and (j^) (which are now written for a single brane located at 
x 5 = 0), one generally finds solutions without imposing any relation among the parameters of 
the Lagrangian. In particular, the brane tension f(4>) - which contains the contribution of the 
fields living on the brane (the SM fields) to the four-dimensional cosmological constant - needs 
not be fine-tuned and may vary over a wide range of values. 

However, these solutions either do not localize gravity on the brane (the graviton zero mode 
is not a normalizable bound state, or equivalently the four-dimensional Planck mass diverges) 
or they have curvature singularities at finite proper distance^. At these singularities, the warp 
factor vanishes - implying that the four-dimensional metric degenerates, — > - and the 
curvature scalar diverges, signaling a breakdown of Einstein gravity. In the absence of any 
mechanism that would smooth these singularities, one assumes that their effect is simply to cut 
off the fifth dimension, effectively compactifying it to a finite interval and ensuring localization 
of gravity on the brane. However, due to the presence of singularities the consistency condition 
( lOl) is not satisfied, which indicates that the ansatz (|Tl|) does not solve the equations of motion 



on the whole interval. This can be seen in explicit examples by the fact that the effective 
four-dimensional cosmological constant, Eq. (|), 

XM 2 Pl = - 2 -Jdx"e' A Vm) - \e iA fm)\ x5=0 (12) 

does not vanish, in spite of the Poincare invariance of the geometry on the brane. Thus, in order 
for the solution to make sense, the (unknown) microphysics that resolves the singularity has to 
contribute to A in such a way that it exactly cancels the brane and bulk contributions. This 
non-trivial statement ruins the apparent self-tuning behaviour of the solution: the traditional 
fine-tuning problem of the cosmological constant has simply been shifted into an unknown 
sector of the model. Of course one cannot exclude that an adjustement mechanism, possibly 
involving some new light degrees of freedom p0|, |16||, may restore the self-tuning; but a concrete 
realization of this idea is still missing. 

In order to make the problem explicit, one can parametrize the contribution of the unknown 
physics at the singularity by adding delta source terms in the energy-momentum tensor, Tmn = 
— J2 S fs(<P)^( x5 — x s )gfj, u 5^f5pf + . . ., where x s denotes the position of a singularity. The necessity 
of adjusting precisely these terms for each particular solution (for each value of the brane 
tension) spoils the self-tuning. Indeed, in addition to the boundary conditions at the brane, 
the solution has to satisfy, at each singularity x s , the following two constraints: 

e iA ^[<f ) '(x s + 0)- ( j ) '(x s -0)} = \z AA{Xs) ^tWs)\, (13) 

e 6A ^ [A'(x s + 0) - A'(x s - 0)] = -\ e 6A ^f s [<P(x s )} . (14) 

o 



3 The authors of Ref. |T(| have shown that in this class of models, the only solutions that localize gravity 
with an infinitely large extra dimension involve a fine-tuning between bulk and brane parameters, like in the 
Randall-Sundrum model. 
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The fine-tuning implied by Eq. and ([14]) guarantees that the consistency condition (|II]) is 
fulfilled and, as a consequence, that the 4d cosmological constant vanishes. 

In Ref. [[L2|], this point was illustrated in the case of a vanishing bulk scalar potential, 
V((f>) = 0, addressed in Ref. and fTOfl . In this section, we would like to discuss another 



example and to show that the requirement of having a global solution of the five-dimensional 
equations of motion in the presence of curvature singularities imposes a fine-tuning for more 
general ansatze than the (four-dimensional) Poincare-invariant background (|TT|). 

Let us first consider the case of a flat brane located at x 5 = 0, i.e. we take the Poincare- 
invariant ansatz (|TT|). In order to illustrate our point, we choose examples of bulk scalar 
potentials for which one can obtain a simple analytical expression for the solution of the equa- 
tions of motion. The simplest case corresponds to a vanishing bulk potential and has been 
addressed in Ref. fll2|. Let us repeat shortly the discussion here. The solution of the bulk 



equations reads [|TU| 



</>(x 5 ) 



[ feilog 


4™5 
3 X 


+ Ci 


+ <*1 , 


x 5 < 


\ ! £2 log 


4-5 
3 X 


+ c 2 


+ d 2 , 


x 5 > 



(15) 



A' (x 5 ) 



Ei (j)' (x 5 ) , X 5 < 

e 2 <p' (x 5 ) , x 5 > 



(16) 



where ei i2 = ±1 and C\, c 2 , d\ and d 2 are integration constants. The boundary conditions 
on the brane determine three of them in terms of the fourth one and of the parameters in 
the Lagrangian; this means in particular that there are solutions for a wide range of values 
of the brane tension and of the scalar coupling to the brane. The x 5 -dependence of the warp 
factor (e 2j4(x ) cx ||x 5 + q| 1//2 ) does not allow for a localization of gravity on the brane with an 
infinite extra dimension; therefore, the only phenomenologically acceptable solutions are the 
ones that have singularities on both axes x 5 < and x 5 > 0, effectively compactifying the 
fifth dimension to a finite interval. In practice this means that we must choose C\ > and 
c 2 < 0; the singularities are then located at x_ = — fci and x + = — |c 2 . Following Ref. fl2]|, we 



truncate the solution at the singularities, i.e. we set |x — x_| =0 for x < x_ and \x — x + \ = 
for x 5 > x+. As explained above, (|T^) and ( |T5D do not solve the equations of motion on the 
closed interval [x_, x + \. This is reflected in the fact that, as can be immediately seen from Eq. 



(|12|), the effective four-dimensional cosmological constant does not vanish, unless the brane 
tension itself vanishes. A simple way to cure this inconsistency is to "resolve" the singularities 
by adding the following source terms to the action (fj): 



- d A x 



-9 f -(<!>) 



- (Fx- 



-9f + (<f>) 



(17) 



The matching conditions (|13D and (0) then amounts to a fine-tuning of the source terms at 
the singularities. Assuming for example exponential couplings for as in || |10[, /(0) = e^T, 
f _((/)) = e fe -*T„ and f + (<j>) = e b +<f>T + , one obtains 



&- = -! 

b + = +l 



-2e + 3 dl 
-2e~l d2 
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and e x = — e 2 = +1 for |6| < | (solution (I) of Ref. HlQj), and 



4 

'3 ' 



T 
' 2 ' 



(19) 



and ei = €2 = +1 for |6| = 
that since in the case \b\ 



[case considered in Ref. |§ and solution (II) of Ref. UlO|). Note 
I the solution is symmetric under x 5 «-> — x 5 , it is possible to 



identify the two singularities and to treat x 5 as a periodic coordinate (this amounts to continue 
periodically the solution beyond the singularities instead of "truncating" it as we did above). 
Then one needs to add a single energy source at x s = x^ = x + , — J d A x \[^g e bs ^T s \ x $ =Xa with 
6, = —-5 and T s = — T. 



The second explicit example |M| we would like to discuss corresponds to a constant bulk 
scalar potential, i.e. a nonvanishing bulk cosmological constant, V((f>) = A#. For the sake 
of simplicity, we focus on the case A# < and consider only Z 2 -symmetric solutions; when 
convenient we shall work on the half line x 5 > 0. The solution of the bulk equations reads [21] 



,4A(0) 



cosh(co>|x 5 |) + r sinh(a;|x 5 |) 



(20) 



(f>'{x b ) 



ce 



-4[A(a; 5 )-A(0)] 



(X 



c I dye 
h 



-4[A(v)-A(0)] 



+ d , 



16 



r 2 - 1 



(21) 
(22) 



where u = J — 3 Ab- Note that the last equation requires \r\ > 1. The boundary conditions on 
the brane give the following two relations 



mm 



-3oor 



df 



[0(0)] 
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(23) 



Plugging (p3|) into (p2[) , one obtains a relation between f [(f>(0)], §^ [0(0)] and A#. Given f(<f>) 
and As, it is generally possible to find a value of d = 0(0) which satisfies this constraint for 
a wide range of values of the brane tension f [(f)(0)]; Eq. (^) then determines the other two 
integration constants r and c. For instance, if has only exponential couplings, f((f>) = e b< ^T, 
(p(0) is determined to be 0(0) = ^ In 16o; 2 /T 2 (y — b 2 ) (note the restriction on the parameter 

b, \b\ < g). Therefore, this solution has a self-tuning behaviour; however, it localizes gravity on 
the brane only if r < —1. In the case r < —1, there is a singularity at finite proper distance, 
located at 

x s = — Arctanh ( — ^ (24) 
(together with its twin singularity at x 5 = —x s ), and the solution can be rewritten 



3 4A0 5 ) 



,4A(0) 



sinh [u(x s — I or 
sinh(co>x s ) 



(25) 
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tanh 


2 \ x 


s 


x 5 \)' 


tanh 
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(26) 
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where e = sign f|^[</>(0)]). In the case r = —1, the solution becomes e iA ^ = e 4A ^e ^' x ' 5 ' and 
the singularity is pushed to infinity. This limit corresponds to a constant <p background (i.e. 
4>(x 5 ) = est, as implied by Eq. (|22|)), and Eq. ( |23| ) amounts to a fine-tuning between the bulk 
cosmological constant A# and the brane tension Vq = /[0(O)], namely Vq = 3u. It is interesting 
to note that this fine-tuning is precisely the one appearing in the Randall-Sundrum model 
(here the configuration would be the one considered in Ref. ||, in which the second brane is 
pushed to infinity; note that our notations correspond to setting 2M = 1 in the formula of 
Randall and Sundrum). This difference of behaviour of the two solutions (fine-tuning versus 
self-tuning) is already striking, given the fact that the fine-tuned solution can be obtained from 
the "self-tuning" solution by taking the limit x s —>■ oo. If one now computes the effective 
four-dimensional cosmological constant, one finds, using fll2]), 

sinn(wi s ) 

Thus, in the case r < — 1 ("self-tuning" solution with a singularity at finite proper distance), 
there is only a partial cancellation between the bulk and the brane contributions to the 4d 
cosmological constant, while in the case r = — 1 (fine-tuned solution with a singularity at 
infinity) sinh(ix>x s ) = oo and therefore A = 0. This resolves the apparent paradox noted above: 
the "self-tuning" solution is incomplete, i.e. an additional source has to be added at the 
singularity in order for the equations of motion to be globally satisfied. Specifically, adding 
— / d^x \/—g e bs ^T s \ x 5 =Xa to the action, one finds that a fine-tuning of b s and T s is required, 

4 m 3 ue~ bsd . . 

bs = r£, T s = - ^ . 28 

3 2 tanh (*f-) 



The contribution of the singularity to the 4d cosmological constant then exactly cancels (p7|). 
as expected. In the limit x s — > oo, this contribution vanishes. 

Another way to understand the necessity of adding an energy source at the singularity is 
to consider a slightly different set-up, in which a second "end of the world" brane is placed at 
a position x 5 = L < x s . Then Eq. (|23| ) must be supplemented with a set of new boundary 
conditions at x 5 = L, 

f L [<l>(L)] = -3cu coth[u{x s - L)] , (29) 

umi 16 sinh(^x s ) 

— - [4>{L)\ = - — c -t • 30 

o<p 3 sinh [u(x s — L)\ 

Since the integration constants d, r (hence x s ) and c are determined by the boundary conditions 
at x 5 = 0, Eq. (|29D and (|30D imply a fine-tuning of /l(0), which at the same time fixes the 
inter-brane distance L, in the same way as in Ref. P5| . More specifically, for an exponential 
coupling of the bulk scalar field f L (<p) = e bL</) T L , Eq. fl29|) and (0) yield 

4 

b L = - e cosh -1 [uj(x s — L)\ , (31) 



( tanh ^(x s — L) 



cosh 1 [uj(x 3 — L)] 



T L = -3ue^ d coth[u;(x s -L)}^ Jw^x L ] (32) 
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Thus, for a given value of bi (resp. T L ), Eq. fl3T|) (resp. Eq. fl32|)) determines the inter-brane 
distance L, while Eq. fl32|) (resp. Eq. (|3T|)) tells us that the brane vacuum energy T L (resp. the 
scalar coupling 6 L ) must be fine-tuned. Conversely, for a fixed value of L, both i>£ and T L must 
be fine-tuned. If one now sends the second brane to the singularity, L — > x s , one finds exactly 
bi — > b s and T L — > T s : the "self-tuning" single-brane set-up is nothing but the (singular) limit 
of a fine-tuned two-brane system, in which the second brane hides at the singularity^. 

Note that we could "screen" the singularities in the same way in the case of a vanishing 
scalar potential discussed previously, i.e. we could place two branes at x 5 = L_ > x_ and 
at x 5 = L + < x + . Unlike what we have just found for a constant scalar potential, we would 
then conclude that the brane tensions and scalar couplings do not depend on the positions 
of the two additional branes, i.e. Eq. ( |i~8| ) (resp. Eq. (|TJ|) in the Z 2 symmetric case) hold 
for any inter-brane distances, and L_, L + are moduli. The statement that one recovers the 
"self-tuning" solution in the limit L_ follows then trivially. 

We could go on discussing other examples with a more general bulk scalar potential; however, 
while the explicit form of the solutions would become more involved, the conclusions would 
remain the same. 

We would now like to show that the necessity of adding energy sources in order to ensure 
the consistency of a solution of Einstein equations with singularities is not a particularity of the 



ansatz (|TT|) . For this purpose, we come back to the case of a vanishing bulk scalar potential, in 



which curved solutions with maximal symmetry (de Sitter or anti-de Sitter) in four dimensions 



can be found [ 20[ . Our ansatz is 

ds 2 = e 2A{x5) g^dx»dx u + (dx b ) 2 (33) 

where g^ u is the metric of a maximally symmetric 3 + 1 dimensional space with a curvature 
constant A (the explicit form of g^ u has been given at the beginning of Section |2|). With this 
ansatz, the solution of the bulk equations of motion can be written in terms of a hypergeometric 
function F(z) = 2 Fi{\, §, §, z) 

F (-^ e 6 ^ 5 )) = q(^ 5 + Cj ) , (34) 

0'(* 5 ) = 7l e" 4 ^ 5 ) , (35) 

where i = 1, 2 refer to x 5 < and x 5 > respectively, ei j2 = ±1, c\, c 2 are integration constants 
and 7i = ei-F^^o/Q if we choose A(0) = 0. The boundary conditions on the brane give the 
following two constraints, 

/[0(O)] = 2 ( V71 2 + 9A + V7! + 9A) , ^ [0(0)] = ~ (72 - 71) • (36) 

Thus, given f(<p), one generally finds a continuous set of solutions corresponding to different 
values of A (including a flat solution A = 0). The authors of Ref. and ^(J have shown that 

4 It is interesting to consider another limit, in which the singularity is pushed to infinity, while the second 
brane remains at a finite proper distance L from the brane at the origin. Then one recovers again the Randall- 
Sundrum model (in the configuration of Ref. Q), with the double fine-tuning Vq = —Vl = 3w, where Vl = 
lim^^oo f L [4>(L)]. 
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for an exponential coupling of the scalar field to the brane, there is a particular value of the 
coupling which does not allow for any curved solution. Indeed, pluging /(</>) = e^T in the 
above constraints (|36|) , one sees that only A = leads to a solution in the case b — ±|, while 
any A is allowed for \b\ < | (and any A < is allowed for |6| > |). Here we are interested in 
curved solutions, so we choose b ^ ±| and consider A 7^ 0. The solutions have singularities on 
both axes if we take c\ > and c 2 < 0, like in the A = case. The vacuum energy is given by 
Eq. (0) with an additional term 12KM 2 PI on the RHS; assuming; first that the singularities do 
not contribute, one finds 

2 



XM 2 Pl = -- ( vVi 2 + 9A + v/ 72 2 + 9A ) + 12AM|, , (37) 



while the Planck mass is computed to be 

1 



m pi ~ ~T 
Pl 9A 



(v/ 7l 2 + 9A + + 9a) - (| 7l | + | 72 |) 



(38 



Since R = 12 A, the four- dimensional Einstein equations require A = 6A. Clearly this is not the 
case here, which indicates an inconsistency of the solution. Adding now delta energy sources 
f±(4>) at the singularities, one obtains two new sets of two boundary conditions, 

e 4 V ± (0) x&=x± = -2| 7l! |, (39) 



x°=x± 



T 7T 7|? , (40) 



where 72 = 7l for a; = x_ and 72 for x = x + . It is now straightforward, using Eq. (p9|), 
to check that the singularities give a contribution | (| 7 i| + | 72 |) to XM Pl , therefore ensuring 
A = 6 A, as required by the 4d Einstein equations. Note that, like in the A = case, the 
boundary conditions at the singularities introduce a fine-tuning which is naively absent if one 
does not insist on the validity of Einstein equations at the singularities. Thus even curved 
solutions require a fine-tuning. 



4 More on nearby curved solutions 

In the previous section we have seen that solutions whith singularities are not consistent unless 
one specifies the physics at the singularity. Our proposal was to interprete the corresponding 
curvature singularities as arising due to the presence of additional source terms. This proposal 
passed the consistency condition ([J). The lesson we learned is that fine tuning on the parameters 
of sources at the singularities is needed. Therefore, these setups do not solve the cosmological 
constant problem. However, irrespective of the fine-tuning problem, they possess an additional 
property which makes them attractive in view of the cosmological constant problem: for the 
very specific choice of parameters (namely vanishing bulk potential and particular value of the 
scalar-brane coupling) discussed in and in solution (II) of [[TO], there are no nearby curved 



solutions U |2(|. In the present section we are going to rederive that result in a simpler way, 
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and afterwards we will explore what happens when one adds a bulk potential for the scalar 
field 0. 

The action considered in |9|, [l(J reads 



S = J d 5 x\^G 



/ d A x^g~Te h 

Jx 5 =0 



(41) 



In our previous notations (Q), this corresponds to V(4>) — and /(</>) = Te b< ^. This is some 
version of five dimensional Jordan-Brans-Dicke (JBD) theory in the Einstein frame where the 
scalar <ft couples in a universal manner (as a conformal factor of the metric) to matter. However, 
all matter is confined to live on a 4d hypersurface located at x 5 = 0. The metric on this 
hypersurface is taken as the induced one. We define the JBD frame by performing a Weyl 
transformation such that decouples from all matter. This is achieved by the field redefinition 
Gmn —> e~ b ^GMN- In the JBD frame the action reads, 



S = J (fxV^Ge 



d x\J—gT. 



(42) 



Now, for b = ±3 we observe that the </> equation of motion results in the constraint R = and 
we see that this case is special. It corresponds exactly to the setup where only flat solutions 
exist || |2D|. Now, we will re-derive this result in the JBD frame (fl"2]) where the calculation 
turns out to simplify substantially. This will allow us to study how the solution is affected by 
the presence of a bulk scalar potential. We should remark that in this section we understand 
that, in the spirit of the previous section, sources at singularities need to be added. However, 
in order to keep formulas short we do not write down those terms explicitly. Typically we are 
looking for solutions of the form (fi,u = 1, . . . , 4) 



ds 2 



^ A ^g^ (x»)da?dx" + 



dx" 



and 



= 0(x 5 ) 



(43) 



(44) 



(a non trivial 55 component of the metric which may appear after the Weyl transformation 
can be absorbed in a coordinate transformation). Further, we chose g^ v to be the metric of a 
maximally symmetric 3+1 dimensional space 



(45) 



where the bar indicates that the Riemann tensor is computed with respect to the metric g^, 
whose explicit form was given in Section ^. 

The equations of motion derived from (|42|) for |6| = | are (as in the previous sections, a 
prime denotes derivative with respect to x 5 ) 



3Ae~ 2A + 2A" + 5A' 2 = 



QAe~ 2A + QA' 2 - ^A'(f>' = 



(46) 
(47) 
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3Ae~ 2A + 3A" + 6A 



l2 



," + 0' 2 + -Te^S (x 5 ) = 0, 



(4£ 



where (|46f) is the equation of motion, and Q47Q , ( f48|) are the 55, and /xz/ components of the 
Einstein equations, respectively. We are looking for solutions with continuous functions A and 
<p. The advantage of the JBD frame is that the source term appears only in one of the equations 
. From we see immediately that A' must be continuous such that A" is finite. With 
follows that either <$' is also continuous or A'(0) =0. If 0' is continuous then fl48|) implies 
that T has to vanish - a case in which we are not interested. Therefore we take A'(0) = which 
restricts the value of A to vanish. So, without solving explicitly any differential equation we 
rederived that for |6| = | there are no nearby curved solutions. Just for completeness we give 
the solution to the remaining equation fl48Q. Solving in the bulk gives 



A(x 5 



const 



U E 

-—log OiiX + Pi 

\ b \ 



(49) 



where the index i = 1, 2 labels the case of x 5 < 0, x 5 > 0, respectively. Using the !Z 2 symmetry 
of the modulus function we can without loss of generality restrict to the case that a« > 0. 
Further, we want to have a finite 4d Planck mass, and this forces us to consider a solution with 
a singularity appearing for some value of x 5 . Then continuity at zero leads to, 



Pi = ~p2 = P > 0. 

The jump condition on the first derivative of (f) at zero finally restricts 



«1 + « 2 



It 

2 



(50) 



(51) 



(for a symmetric solution this implies ot\ = a 2 = jT). Transforming back to the Einstein frame 
and performing an appropriate coordinate transformation in the fifth direction one finds easily 
that this is solution II of flQfl , or equivalently the solution discussed in (i.e. one recovers 



Eq. (H) and (pD with e x = e 2 



-C2 and d\ = d 2 



As a next step we want to investigate how stable the constraint A = is against deformations 
of the bulk Lagrangian. To this end, we will focus on the special case of adding an exponential 
bulk potential for the dilaton 0, insisting, however, on |6| = | . (To some extent this mimics 
also the option of screening the singularities by additional branes - the additional bulk potential 
can be viewed as filling the space with a continuum of three-branes.) To be specific, we modify 
the action (42) as follows, 



S = ( d 5 xV^G le-^R - AeM fe H - / d^x^T 

J l J Jx 5 =0 



(52) 



which corresponds to adding a bulk potential V{4>) = Ae a< ^ in the original Einstein frame. The 
equations of motion with the ansatz (H3f) , (fH) read 



3Ae~ 2A + 2A" + 5A' 2 + 

126 



Ael a -i 







(53) 
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- 6Ae- 2A + 6A' 2 - ^A'(f>' + -Ae( a ~^ = 
\b\ y 2 



3Ae- 2A + 3A" + QA' 2 - 3^A'<j)' - ^0" + f + -AeMV + -Te^8 

\b\ \b\ 2 2 



x 



(54) 



0. (55) 



Together with the requirement of having continuous functions A and the first equation ( |53| ) 
implies that A' must be continuous as well. In order to allow for a non vanishing T in the third 
equation (^) we need a jump in <$' at x 5 = 0. From the second equation (54) one sees that 
this is possible only for 

A'(0) = 0, (56) 



and that any solution has to satisfy the condition 



6Ae -2A(0) = i Ae (a-f>(0)_ 

2 



(57) 



This is quite a remarkable result. Insisting on \b\ = ~ but allowing for non-zero A we find 
that the effective 4d cosmological constant A necessarily differs from zero. The situation is 
somewhat complementary to the case |6| < |. There nearby curved solutions exist even for 
vanishing A [J2DJ. On the other hand it is possible to have flat solutions for non- vanishing A 
(solution /// in and generalizations thereof [16]). To complete the discussion, we now solve 
the equations in the particular case where the warp factor is constant (which trivially satisfies 
the constraint A'(0) = 0), and show that the corresponding solution gives the flat solution (|49| ) 
in the limit A — > 0. Plugging A'(x 5 ) = into the equations (|53|) and (|54"D , one finds that they 
can be solved only if 

a = ~ (58) 

(conversely, if we choose the bulk scalar potential to satisfy (|58|), the equations of motion imply 
that A'(x 5 ) = 0). Condition (|7]) then becomes 



Ae 



-2A(0) 



1a. 
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For A < we find 



— j3 cosh px 5 



ot2 sinh px 
P 



x 5 < 
x 5 > 0, 



(59) 



(60) 



where p 



In order to have singularities we chose ati, (3 > 0. The function is continuous, 



and the jump conditions on the first derivatives leads to (0). For A > we replace in 
the hyperbolic functions by their trigonometric counterparts and p — v 
corresponds to p —>■ 0, and one easily rediscovers (P^). 



f . The limit A ->• 



5 Further constraints on self-tuning models 

Up to now, we have discussed the self-tuning proposal at the level of an effective classical La- 
grangian. We did not specify the relation of these classical models to a realistic quantum field 
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theory. In particular, making such a relation one needs to specify the meaning of the brane 
tension in terms of the perturbative Lagrangian for matter fields, and any other, nonperturba- 
tive in 5d, pieces of the Lagrangian on the brane. Firstly, as we argue later on, there may be 
some primordial nonzero contribution to the tension on any wall, which is the remnant of the 
higher dimensional theory, and can be considered as nondynamical and nonperturbative in five 
dimensions. Such contributions do not even need to be a part of a supersymmetric sigma model 
on the brane, so they are not protected, or set to zero, by brane supersymmetry. Further, there 
is the matter Lagrangian which has perturbative couplings and which is modified by quantum 
corrections on the brane. This Lagrangian may be globally supersymmetric, or may just be the 
one of the usual Standard Model. There are basically two valid points of view on the place of 
the quantum corrections in the effective classical brane tension. One approach is based on the 
observation that there exist (fine-tuned) Poincare invariant solutions for a wide range of cou- 
plings of the bulk scalar to the wall. Then, even if each loop contribution on the brane comes 
with a different dependence on 0, the expectation value of the sum of them up to any arbitrarily 
choosen order of the perturbation theory can be identified with brane tension in our models 
and efficiently screened. This point of view is valid in the present paper up to the section 4, 
and we take in that part the brane tension to be — e b< ^/ M(4) (< V^j'"^ 1 ' > +T P ), where T p is the 
primordial brane tension. However, in such a situation the Poincare invariant solution is not the 
only one, there exist nearby curved solutions which give nonvanishing cosmological constant in 
4d. In a more ambitious approach, one would demand, that there should be no nearby curved 
solutions. This, as shown in |J, requires a very careful choice of the the brane model and its 
coupling to the bulk scalar. There one has to assume that the classical brane tension is just 
the sum of some primordial tension and the vacuum value of the classical Lagrangian for fields 
living on the brane, and watch explicitly the quantum corrections to see whether the required 
form of the model does not get disturbed. To be more specific, we saw in the previous section 
that if we want to achieve certain uniqueness of Poincare invariant solutions, it is important to 
have a very particular value of b, like |6| = 4/3 of the previous section. It has been argued in 
that it is possible to write down a model where quantum corrections respect the value of b once 
it is set at tree-level. The associated stability of the Poincare invariant solution with respect 
to curved deformations is in fact the second part of the original self-tuning proposal. Hence 
it is instructive to take a closer look at the specific brane model assumed in 0. The action 
integral is J d 5 xy/^G(R - |(<90) 2 ) + J d 5 x^/ = gS(x 5 )L b (H;g lxu ) where = g^e^, g^ is the 
induced metric on the brane and \b\ =4/3. The special feature of this model is the coupling of 
the scalar field to the brane in its equation of motion 



we were in four dimensions, the equation of motion for the scalar which couples to the trace of 
the energy-momentum tensor only would automatically imply vanishing of that trace on-shell, 
hence the vanishing of the cosmological constant. However, in the present case the vacuum 
value of derivatives of the scalar with respect to the coordinate transverse to the brane can 




(61) 
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be nonzero without violating the Poincare invariance. Hence, the equation of motion does not 
demand the vanishing of the trace of the brane energy-momentum tensor. In fact this is the 
situation we consider in this paper. However, the nonvanishing of 9^ implies the scale anomaly 
on the brane which in turn implies that the Weyl rescaling on the brane is anomalous. This 
means that the conjectured decoupling of the scalar from the brane, hence 'conservation' of b 
by brane physics, may be violated in the full quantum brane theory. In this respect, we agree 



with the argumentation put forward in Ref. |JL4||. On the other hand, one could imagine choos 



ing the brane sector which is conformally invariant to all orders, like N = 4 Super- Yang-Mills 
models. This implies 9£ = at all orders and the background cosmological constant vanishes 
as a result of a symmetry. In a realistic model however, such a symmetry has to be broken, 
and the cosmological constant problem reappears. 

It is interesting to write down explicitly the 1-loop corrections to the effective potential 
in an example of a model conformally coupled to the scalar. To this end we recall that the 
ultimate contribution to the vacuum energy which must be cancelled in order to be consistent 
with observation, arises in low-energy theory, say around the 1 TeV scale. At that scale in 
all realistic models 4d brane symmetries which may have something to do with protecting the 
vacuum energy, like N = 1 supersymmetry or conformal invariance, are broken. This is the 
place where the mechanism of self-tuning (if not fine-tuned) would be really useful. To see 
what happens when quantum effects are taken into account, let us specify the representative 
content of the brane to be given as 

e H I v^(e~ 6/2 V^d^ - m 2 ^ + *e"^ (62) 
J Ag 

To compute one-loop contribution to the effective potential we first perform the redefinition of 
the fields so that their kinetic terms are canonical with respect to Minkowski background on the 
brane, and then use the standard formulae with a physical UV cut-off A which we understand 
to be close to 1 TeV . After going back to the original frame one obtains — J y/— g(e b ^(Vo + 
T p ) +Strl (9(A 4 ) + g2^e fe//2< ^A 2 (m 2 — M 2 )) + ..., where one can see terms which are multiplied by 
different powers of e h ^ . Of course, this situation is also manifest in models without conformal 
coupling of the brane to the scalar field. Hence, in general, even if one carefully adjusts 
|fe| = 4/3, at tree level, one departs fom this choice after including quantum corrections. This 
simply means that the classical selection rules have limited value in the true, quantum, world. 

We want to mention another troublesome feature of the model with conformal coupling 
of the scalar to the branes, which is an inconsistency with long range properties of 4d grav- 
ity at the classical level. This becomes clear in the Jordan-Brans-Dicke frame introduced 
in the section 4. In that frame the model of takes the form of a generic Jordan-Brans- 
Dicke model, and is subject to constraints given by the analysis of 4d Jordan-Brans-Dicke 
systems performed in P8| . The 4d Jordan-Brans-Dicke action in the JBD frame looks as fol- 
lows: / cftxyf —g{^R — 2uo{^)g^ h 'd^d u § — V($)) with the assumption that matter fields are 
coupled only to the metric g^ v and not to $. This is very similar to the Lagrangian given in 
(f42| ), up to a redefinition of the field <fi and integration over the fifth dimension. It is easy 
to see, through the integration of the fifth dimension directly in the 5d JBD frame, that in 



the case of (S3 ) uo = 0. It has been shown in 281 that there exist constraints on the value of 
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the parameter u of JBD models. This constraint says that < 10~ 3 , hence in the model 

which we are discussing here the above bound is violated. The troublesome constraint could be 
avoided if the 4d JBD scalar would receive a mass of the order of 1/mm. This can be achieved 
by introducing a very soft potential for the scalar field, so that the experimental bound on the 
4d cosmological constant, p\ < (0.002 eV) 4 is not violated. At this point however, we should 
recall the result of the preceding section, where we have shown that introducing a simple bulk 
potential while insisting on \b\ = 4/3 has led one from a flat to a curved solution. 

However, one should repeat, that in the context of generic models lacking conformal cou- 
pling of branes to a bulk scalar, it is still likely that even for a modified, so pretty generic, 
dependence on there exists a Poincare invariant vacuum background. The problematic point 
is, that in response to quantum corrections the tensions on the additional branes located at 
singularities must get re-finetuned. This would have to be achieved by some mechanism oper- 
ating at low energies in addition to the bulk scalar mediation. One may speculate that the bulk 
KK modes could be an agent of such an readjustment, but a separate calculation is needed to 
decide whether this is the case. 

The next issue which needs to be addressed is the physics at or near branes located at 
singularities. One way to conceal this problem is to say that the singularities are resolved 
in the microscopic higher dimensional, D > 5, model. However, perhaps one can still make 
some useful observations looking more carefully at the 5d models. A way to test physics 
near singularities is to perturb the model by a small point-like test mass and to read off the 
strength of interaction of this test mass with the gravitational field. The action for such a test 
particle is S p = —m J ds in 5d coordinates and it gives a contribution to the energy momentum 
tensor of the form T MN = -7= / dr , x = S(x — x(t)). The perturbation of the action 

under a small variation 5h^ v is 5S P = — J d 5 x^—g(e 2A 5h^ u ) — — J d 4 x^ ' —gbh iiV T iW where 
T^ u = J dye 6A T^ u . One obtains then the 4d energy momentum tensor due to a test mass 

rp^u _ me r ^ x x = S(x — x(t)). This means that the mass which interacts with the 

J ^-gpQxPxQ v v 11 

gravitational field at the location x 5 is me A ( x \ hence, for instance, it vanishes whenever e A = 0. 
Hence, physical matter located at singular branes is effectively massless. Essentially, since 
\J — G = at the singular brane, local dynamics there is also suppressed. Then the mechanism 
which creates correlations between vacuum tensions on physical and singular branes might be 
the result of some higher dimensional consistency condition, rather than a dynamical one. For 
instance, one can recall the way the Bianchi identity is fullfilled in the Horava-Witten model. 
There one needs to switch on vacuum fluxes of the gauge fields along the internal dimensions. 
The fluxes are different on different walls, but their sum is fixed by the Bianchi identity. These 
fluxes, together with that of gravitational curvature tensor in six compact dimensions, lead to 
nonzero but correlated brane tensions in five space-time dimensions. 

To learn something about dynamics in singular warped universe, it is useful to compute 
in the nonrelativistic limit the force acting on a test mass m. From the geodesic equation 
one easily obtains the acceleration of the freely falling particle = —^00 = ~\d MN ^w- 

This is nonvanishing along the fifth direction and gives = A' {x 5 )e 2A ^ x °\ Hence, the 
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force acting on the slowly moving test particle is F 5 = mA'(x 5 )e 2A ( x \ To see the possible 
implications let us take an example of the singular background which has been discussed in the 
third section of this paper, in the vicinity of the naked singularity on the positive side of the 
x 5 axis, e A = e 1 ^ 3d \A/3y — c] 1 ^. When the particle approaches the singularity from the left, it 
feels the force F 5 — |4/3y — c| -1 / 2 which repells the test particle away from the singularity. 
This may be interpreted as a hint at sort of stability of the system, in the sense that one does 
not expect the flow of physical matter from the positive-tension brane to the singular negative 
tension branes. 

Finally, let us comment briefly on the stabilization of the fifth dimension in the static (flat) 
case. The standard way to formulate this problem is to separate a ^-independent part of 
G55, let us call it here r c (x), and to require that after integrating out all remaining degrees of 
freedom and integrating over x 5 there remains an effective potential for r c such that either it 
has a minimum (and thus truly stabilizes r c ) or it is flat, so no effective force moves r c around. 
This second possibility would correspond for instance to a modulus in a supersymmetric model, 
but we know that at low energies supersymmetry needs to be completly broken down. Hence, 
of actual interest is the first one. The examples of that favourable situation are the models 
(with repaired singularities) presented in this paper which contain a nontrivial bulk potential 
for the scalar field. In the case of the models with no potential in the bulk, the distance to the 
singularities was not fixed. Moreover, if one puts the brane screening the singularity between 
the SM brane and singularity, its position is also undetermined. According to what has been 
said earlier, we expect quantum corrections to modify the classical solution in these cases. 
Parenthetically, we would like to notice that in the 5d Einstein frame in which we work in the 
initial sections of the paper, the substitute of the equation of motion determining r c is simply 
the Einstein equation E 55 — |T 55 = 0. In addition, for the general class of solutions discussed 



here, it has been shown in Ref. |27| that they are stable in the sense that there are no tachyonic 
fluctuations. 

The remaining constraint on singular warped universes which we want to mention arises 
when one computes the corrections to Newton's law due to the exchange of heavy gravitons 



2fJ. When one tries to send the singularities far apart, i.e. when one considers a macroscopic 
fifth dimension, then the correction to the effective 4d Newton's law is 5V = — -\fj,Gp{ ^ C1 ^ C2 ^ ■ 
Hence, the correction grows linearly with the size of the fifth dimension, similarly to what would 
happen in a toroidal compactification. This means that the distance between singularities 
should be smaller than, say, 1 mm. In this regime the correction due to heavy gravitons is 

7rr I I 

exponentially suppressed as 5V = — ^/iGatc 2 d c ii+i c 2i) (1 + sin( ,^V , )). 



6 Summary and Conclusions 

In the present paper we have argued that in brane worlds the cosmological constant problem 
is merely seen from a new perspective. The degree of fine-tuning needed to obtain a value in 
agreement with observational bounds is not improved. We gave explicit consistency conditions 
on warped compactifications. These consistency conditions were derived from the requirement 
of having a globally well defined solution to the equations of motion. If the conditions are 
violated the metric on the brane is not at a stationary point of the 4d effective action. Various 



17 



examples with singularities are inconsistent unless the singularities are screened by additional 
sources. The vacuum energy at these sources needs to be fine-tuned. 

Another issue appearing in the literature is that for |6| = 4/3 and vanishing bulk potential 
only solutions with zero cosmological constant exist |], ^(J . We pointed out that this advantage 
turns into a disadvantage as soon as a bulk potential for the scalar is switched on. Then 
only solutions with a non-zero cosmological constant exist (the numerical value is given by 
parameters of the bulk potential). Hence, in addition of fine-tuning \b\ = 4/3 one needs to 
fine-tune parameters in the bulk potential even at the classical level. 

Finally, we pointed out that quantum corrections due to the theory living on a brane may 
alter the functional form of the coupling of the scalar to the brane. Other open problems are 
the stabilization of the scalar and the inter-brane distance. 
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